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Abstract 

We study effects of CP violation in the associated production of a charged 
Higgs boson and a top quark at the LHC, pp — ► tH^ + X. We calculate the 
CP violating asymmetry between the total cross section for H + and H~ 
production at next-to-leading order in the minimal supersymmetric standard 
model (MSSM), and perform a detailed numerical analysis. In the production 
the asymmetry is of the order of 20%. The asymmetry in the production and 
any subsequent decay of an on-shell charged Higgs boson is to a good approx- 
imation the sum of the asymmetry in the production and the asymmetry in 
the decay. We consider subsequent decays of to tb, v t t^- and Wh°. In 
the case with — ► tb decay, mainly due to CP violating box graphs with 
gluino, the asymmetry can go up to ~ 12%. 
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1 Introduction 

If a charged if ± boson is discovered at LHC or at any future collider, it would be a 
clear signal for Physics beyond the Standard Model (SM). The next question would 
be which Physics beyond the SM it is - almost all extensions of the SM contain a 
larger Higgs sector and inevitably predict the existence of a charged Higgs boson. 
CP violation (CPV) is a possible tool to disentangle the different charged Higgs 
bosons. The phenomena of CPV is important also because it is believed that this 
is the key to our understanding the observed abundance of matter over antimatter. 
Most extensions of the SM contain possible new sources of CPV through additional 
CPV phases. 

In this note we study CPV in the Minimal Supersymmetric Standard Model 
(MSSM) with complex couplings, being one of the most promising candidates for 
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an extension of the SM. In MSSM the additional sources of CPV are the phases of 
the higgsino mass parameter /i = |/z|e l< ^ in the superpotential, of the gaugino mass 
parameters = \M i \e % ^' i , i = 1,2,3 and of the trilinear couplings Af = \Af\e l< ^ f 
(corresponding to a fermion /) [1], respectively. ( Usually M2 is made real by 
redefining the fields.) From the point of view of baryogenesis, one might hope that 
these phases are large [2]. Although the experimental upper bounds on the electron 
and neutron electric dipole moments [3] constrain the phase of fx, M < O(10~ 2 ) [4], 
for a typical supersymmetry mass scale of the order of a few hundred GeV, the 
phases of the other parameters mentioned above are practically unconstrained. The 
CPV effects that might arise from the trilinear couplings of the first generation A u ^, e 
are relatively small as they are proportional to m u ^ e . The same argument holds for 
the second generation. Nevertheless, the trilinear couplings of the third generation 
At,b,r can l e& d to significant CPV effects [5, 6], especially in top quark physics [7]. 

Recently we studied the effects of CPV in the three possible decay modes of the 
MSSM's charged Higgs boson into ordinary particles [8, 9, 10, 11, 12] H ± — > tb, 
— > z/r ± and H ± — > W ± h°, where h° is the lightest neutral Higgs boson. Loop 
corrections induced by the MSSM Lagrangian with complex couplings lead to non 
zero decay rate asymmetries between the partial decay widths of H + and H~ , which 
is a clear signal of CPV. 

Studying the effects of CPV in the decay H ± — > tb, we found that these effects 
can be rather large and reach up to 25% [8] . This is mainly due to the contribution 
of the loop diagrams with stops and sbottoms, whose couplings are enhanced by 
the large top quark mass. This motivated our interest in studying CPV also in 
the production of H ± at LHC (considered previously in [13] and [14]), where the 
dominant production process is the associated production pp — > H^t + X, which 
proceeds at parton level through the reaction bg — > iJ ± t [22]. This process contains 
the same H^tb vertex and corresponding loop diagrams as the decay H ± — > tb, 
and one would expect that the CPV effects might be of the same magnitude. In 
addition, in the production process there are box graphs that are of the same order. 
These contain additional sources of CPV and must also be taken into account. 

We assume that the charged Higgs is produced on mass shell and we consider 
the production and decay processes separately. In this paper we first study CPV 
in if ± production at the LHC, pp — > H ± t + X, through bottom-gluon fusion in the 
framework of the MSSM, with running top an bottom Yukawa couplings. Then we 
study the CPV asymmetry in the combined process of H ± production and decay 
into tb and z/r ± , with CPV in both production and decay. We present a detailed 
numerical study for the CPV asymmetry induced by vertex, selfenergy and box 
corrections in the MSSM. 

The paper is organized as follows. In the next section we study the subprocess 
bg — > H^t including vertex and selfenergy loop corrections and obtain analytical 
expressions for the cross section and the CP-asymmetry at parton level. In Section 
3 we add the parton distribution functions (PDF's) and obtain the CPV asymmetry 
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of the pp — > /J ± t + X production process. In Section 4 we obtain the asymmetry 
in the case of charged Higgs boson production and subsequent decay. Section 5 
contains the numerical analysis in the MSSM. We end up with a Conclusion and 3 
Appendices, which contain some detailed formulas needed in the analysis. 



2 The subprocess bg — ► tH^ 

We study the following processes connected by charge conjugation 

br(Pb) + g^{p 9 ) — ts(Pt) + H-{p H -) , (1) 

Kipi) + 9 a M — ► *.(Pf) + H + (p H+ ) , (2) 

where r, s and a are colour indices, r,s = 1,2, 3; a = 1,...,8. In the kinematics 
of the processeses we neglect the bottom mass m&, working in the approximation 
ml jm\ ~ ml/m^y ~ 0. However, we keep non zero in the Yukawa couplings, 
where it is multiplied by tan (5 or cot (5. 

The tree-level process (1) contains two graphs - with exchange of a bottom 
quark (s-channel) and with exchange of a top quark (t-channel), see Fig. 1. The 
Mandelstam variables are 

s = (p b + Pg) 2 , i = (p t - p g f = (p b - p H -) 2 . (3) 




Figure 1: The tree-level graphs of the bg — > tH process. 

At tree-level there is no difference between the cross sections of the considered 
processes (1) and (2). An asymmetry due to CP non-conservation appears at one- 
loop level. There are three types of MSSM loop corrections to both s- and t-channels 
that lead to CPV - corrections to the H^tb-vertex, selfenergy corrections on the in- 
line and box- type corrections, see Fig. 2. The first two types, vertex and selfenergy 
corrections, are analogous to those in the decay H ± — > tb. The CPV effects in this 
decay were studied in [8] . Our analysis in [8] showed that the main contribution to 
the CPV asymmetry is due to the vertex diagram with a gluino and the H ± — W ± 
selfenergy graph with a tb loop. These contributions are enhanced by the large top 
quark mass and the colour factor of 3. The contribution of the rest of the graphs 
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with supersymmetric particles is negligible. We will present analytical expressions 
for the vertex correction with tbg and selfenergy correction with tb in the loops 
in the production process. The expressions for the box-diagram contributions are 
rather lengthy and we do not present them analytically, but in the numerical section 
they are taken into account. 

2.1 s-channel amplitude 

The matrix element of the graph with bottom exchange on Fig. 1 (the s-channel), 
including the vertex correction with a gluino in the H~tb-vertex and the H~ — W~ 
selfenergy graph with a tb loop (see Fig. 2) reads 1 

M s = i 9 ju s (p t )\[{y t + 5Y t s )P L + (y b + 5Y b s )P R ]^ b +j g ) + 

+S[f S R 2 R P L + (ft'l - fL L )P R ]}T: r Yu r (p b )e«(p g ) , (4) 

where yt and yt, are the (real) tree-level couplings and the other terms are induced 
by the loop corrections. The principal difference between the bg — > tH ± production 
and the — > tb decay, considered previously in [8], is in the vertex corrections. 
In the production process one of the quarks in the if ± t6-vertex is always off-shell 

- this is the b quark in the s-channel and the t quark in the t-channel. In the 
H ± — > tb decay all particles are on mass shell. This leads to a different structure of 
the matrix elements. The one-loop form factors of the decay repeat the structure of 
the tree-level couplings, whereas in the production there are new terms in addition 

- these terms appear in the second lines of eq. (4) for the s-channel, and eq. (12) 



for the t-channel. For the one-loop form factors in (4) we obtain 

8Y t s = m- g f/ L + tthUll + III) , SY b s = m- g ff R + mtf& , (5) 

fin = i£hR ct , f LR = n l Lm nt n e^ {G A ) mn , (6) 

fH = ~3^~f' LL > III = ^L^L(^)mn > ( 8 ) 

ftm = ~3^~f RR C i ' f RR = ^RrrJ^Rn( G drnn , (9) 
1 In this section details are given on the 7J~-production only. 
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Figure 2: The main sources of CP violation in bg — > if ± t production. 

/ll = -^%(5o + 25 1 ), (10) 

where a w = g 2 / (Air) and the arguments of the Passarino-Veltman (PV) integrals 
are 

c t = C i( m h m H+i § i m h m h m l)i i = 0, 1,2 

^ = 5,(m^ + ,m?,m?), j = 0, 1. (11) 

The full expressions for the mixing matrices, the couplings, as well as the definitions 
of the PV integrals are given in the Appendices A, B and C, respectively. 



2.2 ^-channel amplitude 

The matrix element of the graph with a top exchange (t-channel) in Fig. 1, including 
the vertex correction with a gluino in the H~tb- vertex and the H~ — W~ selfenergy 
graph with a tb loop (Fig. 2), reads 
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M* = ij^ 2 u s (p t )e«(p g )rT: r Uj> t -j> g + m^ 

+(t - m 2 t )[{ff L + ~f LL )P L + ff R P R ]\u r {p b ) , (12) 

where the one-loop form factors are analogous to those of the s-channel 

8Y} = m~ g ff L + mtifil + III) , 8Y h l = m~ g fl R + m t ff R (13) 

f'A = f B AiCl - Cl), etc., (14) 
but with different arguments in the PV integrals 

Cl = Ci(i,m% + ,ml,m 2 § ,mlml), i = 0,1,2. (15) 

2.3 Cross section - parton level 

In general, the differential cross sections for the processes (1) and (2) are given by 

where rfcr ± are averaged (and summed) over initial (and final) colour and spin of 
all particles of the process, we have used S^ r=1 £® =1 T S "(T S "*) = 4, and the ± signs 
stand for H ± production. 

We write the matrix elements of the processes (1) and (2), including the vertex 
and selfenergy corrections, in the form 

M ± = M tree,± + M loop,± . (17) 

Here J\4 tree < ± = + are the tree-level matrix elements (proportional to y t 
and y b ), and M loop ' ± = M 3 ^ + + M 3 ,^ + are the loop contributions. 

s(t) 

Aii have the same structure as the tree- level matrix elements and are proportional 
to 5Y t s jf. Ai & 2^ contain the additional terms in (4) and (12). For the squared matrix 
elements lA^I 2 , up to terms linear in a s and a w , we obtain 



lA^I 2 = \M tree ^\ 2 + 2Re^( y M tree ' ± )*M locrp ' ± ^ . 



Further, we need to sum over the polarizations of the incoming gluon. At loop 
level special care must be taken to preserve gauge invariance. We use the axial 
gauge: 

a) = ^( - ^ - ^ + ^-y^ ) , (19) 
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where r] is an arbitrary four- vector that fixes the gauge and fulfills r/.e = and 
rj.k 7^ 0. One can see that in this gauge the unphysical longitudinal degree of free- 
dom manifests itself by an ^-dependent polarization sum over the two transverse 
gluon polarizations. The cross section, being a measurable quantity should be gauge 
invariant and therefore the ^-dependence should ultimately cancel. Using only — g^ v 
on the right side of eq. (19) is sufficient at tree-level, because the term resulting 
from the second and third term on the right side of eq. (19) drops out in the squared 
matrix element of the sum of the s- and t-channels. This is not true anymore at 
one-loop level because the one-loop factors, e.g. SY t s b , eq. (5), and SY^ b , eq. (13), 
are different. 

For the terms in (18) we obtain 
M S *'~M S {- = -32na s (5Y t s y t + SYM^s, t) - 2U(s, t)c v ] , (20) 

Ml*'~ Ml'' = -327ra s (5Y t s y t + 5Y b s y b ) [X 12 (s, t) + 2U(s, t)c v + U(s, t)V(s, t)] , (21) 

Mg'-Mt = -32Tra 8 (5Y?y t + SYfa) [X 12 (s, t) + 2U(s, t)c v + U(s, t)V(s, t)] , (22) 

M*'-M\- = -32na 8 (5Yfy t + SY^y b )[X 2 {s, t) - 2U(s, t)c v - 2U(s, t)V(s, t)] , (23) 

Ml*'' MT = 327ra s m t [(ff L - f LL )y t + f/ R y b ] [1 + 2c,] , (24) 

Ml*'' Mr = 32na s m t \{ff L - f LL )y t + f s RR y b }[y(s,t) - 2c v - V(s,t)\ , (25) 

Ml*'~MT = 327ra s m t [(/£i + f LL )y t + ff R y b ][l + 2c,] , (26) 

Ml*'-Mr = 32na s m t \{f t r L + f^Vt + /^P(M) - 2c, - V(s,t)\ . (27) 
Here X, U, V and y~ functions are given by 

X 1 (s,i) = t -^t, (28) 
s 

X 12 (s, i) = (i-^is + ^- rnU-rnts 

s(t - mf) 

Ms,t)- j—^ 2 , (30) 

m(m*) = 5+£ : <+ , v(8,i) = yjh^4, (si) 



s t — m 



y(s,i) = s + m »\ * . (32) 

t — mf 



The terms proportional to the number c v carry the 77-dependence. The calculation 
is done for two specific choices of r] = p b ; c v = and i] = p b + p g ; c v = 1. In the 
center of mass system, p b + p g = 0, it is easy to see that in both cases the above 
conditions 77. e = and rj.k 7^ are fulfilled. The sum of the products (20) - (27) is 
independent on 77 and therefore gauge invariant. 

Using eqs. (16) and (18) for the total parton level cross sections of the processes 
(1) and (2) we obtain 

«r± = a tree + 7-^7^7 f maX Re{(M tree ' ± )*M loop ' ± ) dt , (33) 
8ns 2 96 Jt min I J 

where the integration limits are given by 

t min ,max = 7j ( m t + m H+ "S=F A V2 (S, m%, 771^+)) , (34) 

and a tree is the tree-level cross section, which is the same for (1) and (2) [15] 
° tre Ts) = ^(^ + ^{A 1/2 (s,m^m^0(7(m^ + -m 2 )-3s)- 
-2^2(m 2 H+ - m 2 ) 2 - 2(m 2 H+ - m 2 )s + s 2 ^j x 



s — (m 2 TT+ — m 2 



- \ 1 / 2 (s,m 2 ,m 2 H+ ] 



X ^ I s — (m 2 H+ — m 2 ) + A 1 / 2 (s, m 2 , m 2 H+ ) ) }' 

We write the cross sections of the conjugate processes ^ given with (33) as a 
sum of CP invariant and CP violating parts 

a ± = a inv ± a CP , (36) 

where the CPV part a CP is given by 

& = (X 1 + X 12 +UV)dt 

+ / A\X 12 + *2~ UV) dt- (B s + B t )(l + y-V)dt\, (37) 

with 

A"' 1 = Im (f RL )y t + Im (f LR )y b ] Im (cf>) + ^[ Im (f LL )y, + 

+Im (/„»)» ]Im(cf») -^lM h ,)y, lm(Bo t 2Bl) ,m 

on m w 

B s(t) = ^[ Im (h L )yt + Im (/ M )y 6 ] Im (C$>) , (39) 

and the CP conserving part <7 irw can be expressed in terms of a with exchanging 
the imaginary parts of the couplings and the PV integrals with real ones 

Qinv = & tree _ & CP j m _^ Rg ) _ ^ 
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2.4 CP violating asymmetry - parton level 

We define the CPV asymmetry at parton level as the difference between the total 
number of produced H + and H~ in bottom-gluon fusion 

A cp = -> tH+) - ajbg -> tH-) 
p a(bg -> tH+) + <7(fy -> tif") " { ' 

Taking into account (36), we obtain 

frCP h cp 

A CP = - ~ (42) 

where a CP and a tree are given by (37) and (35). 

3 The LHC process pp — ► tH± + X 
3.1 Cross section 

We study charged Higgs boson production associated with top quark production in 
proton-proton collisions 

p (P A ) + P (Pb) - t(pt) + ^{pni) + X. (43) 

The Mandelstam variable is S = (P A + P B f ( for LHC VS = 14 TeV). We set 
Pb = x\)Pa = ibPp and p g = x 9 Pb = x 9 Pa, where Xi (x~i) is the momentum fraction 
of the hadron B(A) carried by the parton i. Neglecting the proton mass compared 
to \/~S we get s = XbX g S = XbX g S. We have 

a~(pp — > tH~) =2 fb(xb) f g (xg)dr'(xbXgS)9(x b XgS - S )dx b dx g , (44) 

(J+(PP _> tH + ) = 2 j f~ b (X- b ) ^ fg(Xg)d- + (X- b XgS)6(X- b XgS ~ S q) (IX^X g . (45) 

Here So = (m t + m H +) 2 fixes the kinematically allowed energy range, and /& and 
f g are the PDF's of the bottom and the gluon in the proton. As fb(xb) = fi(xi)i 
we obtain 

a ± (pp -> tH~) = 2 J fb(x b ) fg{xg)a ± (xbXgS)9(x b XgS - S )dx b dx g . (46) 

The factor 2 in the above expressions counts the two possibilities -b (g) comes from 
the proton A (B) and vice versa. 
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3.2 CP violating asymmetry 

We define the CPV asymmetry at hadron level as the difference between the total 
number of produced H + and H~ in proton-proton collisions 

a cp = viPP -> { H+ ) ~ viPP -> tH ~) (A7) 
p a(pp -> tH+) + a(pp -> tif-) ' 1 J 

Taking into account (46) we obtain 

a cp = [ fb(xb)fg(xg)(a + - g-)9(x b x g S - S )dx b dx g 
P I fb(x b )fg(xg)(a+ + a-)6(x b x g S - S )dx b dx g ' 

According to (36), for the CPV asymmetry Ap P up to terms linear in a s and 
a w , we obtain 



CP 

A C /=°-^, (49) 



(jtree 



where a CP is the CPV part of the cross section 



a CP [pp -> tH~) = 2 / / 6 (x 6 ) / f g (x g )a CP (x b x g S)e(x b x g S - S )dx b dx g , (50) 
Jo Jo 

and <r tree is the tree-level cross section 

a tree (pp - fif") = 2 jf 1 / b (z 6 ) jf 1 ^(^^^(x^^)^^^ - S )dx b dx g . (51) 



4 H^ 1 production and decay at LHC 

After the charged Higgs is produced in proton-proton collisions it will be identified 
through some of its decay modes. Here we study the combined processes of H^ 1 
production and decay considering H ± decays into tb, and W ± h°. 



4.1 The subprocess bg -> tE ± -> tt'b (tu^; tW ± h°) 

In the narrow width approximation, when the decay width of if ± is much smaller 
than its mass m H +, the total cross section for charged Higgs production in bg — > 
tH ± , with a subsequent decay — > /, where / stands for the chosen decay mode 
/ = tb; vr^- and W ± h°, is given by 

<r± = <7 P (&5 - tH±) -±= J X (52) 

1 H+ 

Here, <7p is the total production cross section, r(if ± — > /) = is the correspond- 
ing partial decay width of H ± 1 and Y H + is its total decay width. 



11 



We already had the expression for the production parton level cross section a$ 
in the form (eq. (36)) 



(7, 



±a 



CP 



(53) 



The considered partial decay widths of H ± , assuming CPV, were obtained in [8, 
9, 10], and we write them in the form 

rf = r f (rf v ± r^ p ) , (54) 

where V 1 ™ and Vj p are their CP invariant and CP violating parts. 

For the total cross section of if ^production and decay at parton level, assuming 
CPV in both production and decay, we obtain 



1 



Up 1 j- ± 



a P 1 j- +1 



/ 



(55) 



4.2 CP violating asymmetry - production and decay 

We define the CPV asymmetry in charged Higgs boson production in pp — > tH ± : 
with a subsequent decay H ± — > /, assuming CPV in both production and decay, as 

a cp = <PP ~+ tH+ -> if) - a{pp -> tH~ -> t/) 
/ a ( pp ^tH+^if) + a{pp^tH-^tf)' 1 J 

where / stands for the chosen decay mode / = tb; ut + and W + h°. 
In narrow width approximation, taking into account (52), we get 

a cp = *P(PP ~+ iH+) Y{H + -> /) - <r P (pp -> t/f) r(i7~ -> /) 
/ ap ( pp ^tH+)T{H+^f)+op{pp^tH-)T{H-^f)' 1 j 

which leads to [13, 14] 

^inv-pCP I -CPri™ ^.CP pCP 

/ -.invT^inv ^inv ~pinv P D,f ' V / 

Cp i j ^ P f 

i.e. when the decay width of H^ 1 is much smaller than its mass m H +, the total 
asymmetry A^ p is an algebraic sum of the CPV asymmetry Ap P in the production, 
and the CPV asymmetry Ap, P f in the decay / of the charged Higgs boson. 2 



2 In [16] the asymmetry A^ p in the decay is denoted with S CP . 
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5 Numerical analysis 



We present numerical results for the charged Higgs rate asymmetries Ap P , Af b p 
and A„f, eqs. (47) and (56), in the MSSM. All formulas used in the numerical 
code are calculated analytically, except for the box contributions, which are rather 
lengthy. Furthermore, all individual one-loop contributions are checked numerically 
using the packages FEYNARTS and FORMCALC [17]. We also use LOOPTOOLS, 
see again [17], and FF [18]. In the numerical code the Yukawa couplings of the 
third generation quarks (h t , hf,) are taken to be running [8], at the scale Q = 
rriH+ + m t . For the evaluation of the PDF's of the bottom quark and the gluon, f b 
and f g , we use CTEQ6L [19], with leading order PDF's and next-to-leading order 
a s , at the same scale Q. We assume the grand unified theory relation between 
Mi and M 2 , so that the phase of Mi = 0. Our numerical study shows that the 
contributions of the loop diagrams with chargino, neutralino, stau and sneutrino to 
the considered CPV asymmetries are negligible, and besides one exception we show 
only the contributions from diagrams with tb and g. We start from the following 
reference scenario: 

tan/? = 5 , M 2 = 300 GeV , m g = 727 GeV , M = M Q = M D = 350 GeV , 
H = -700 GeV , \A t \ = \A b \ = 700 GeV , <I>m = \, <t>A b = ^ = 4>z = ■ (59) 

The relevant masses of the sparticles for this choice of parameters, and also for 
tan (3 = 30 are shown in Table 1. 

As we will see, in such a scenario the effects of CPV are substantial. It is not one 
of the commonly used minimal supergravity or constrained MSSM scenarios [20], 
for which most studies have been done. There exist experimental constraints from 
b — > S7, relic density, etc.. In principle, there are enough free parameters in the 
general complex MSSM to be compatible with all data. 
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Table 1: Masses of the sparticles (in GeV) for the parameter set (59). 

In Fig. 3 the tree-level cross section a(pp — > tH~ + X) is shown as a function 
of 171h+, based only on the parton process gb — > tH~ , with on-shell and with 
running h t and ht,- Taking running h t and hb reduces a tree by about ~ 30%. For 
m H + > 1000 GeV the cross section drops below 1 fb. 
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a(pp ->■ tH + X) 




400 600 800 1000 1200 



m H+ [GeV] 



Figure 3: The hadron tree- level cross section for the chosen set of parameters (59) 
as a function of m#+ . The green dashed line corresponds to the result using on-shell 
and the solid blue line to that using running h t and h^. 

5.1 Production asymmetry 

As expected, the CPV asymmetry in the production due to the loop corrections 
with tb and g is of the same order of magnitude as in the case of the decay — > 
tb [8], and can go up to ~ 20%. Moreover, the contributions of the box graphs are 
significant and can be dominant for relatively small m H +. 

In Fig. 4 the dependence of the different contributions on of the underlying 
parton process gb — > tH~ is shown for m#+ = 700 GeV. The kinematical threshold is 
at \/I = 871.4 GeV. First the box contribution is the biggest one with a maximum 
at ~ —23 %. Then it drops down asymptotically to zero at \/I > 1700 GeV. 
The vertex contribution has a similar shape with about half of the size of the box 
contribution. But for \f§ > 1500 GeV it becomes constant being of —2%. The 
selfenergy contribution is independent of v§, about —12 %. We show it only for 
completeness. Here also the box and vertex contribution with denoted by 

the pink dash-dotted line, are shown. This contribution is always below 0.5%, and 
therefore we will not show it anymore in the figures. The two spikes in the box and 
in the vertex contributions denote the two thresholds \/I = rrig + mj )i2 , see Table 1. 

The contributions of the vertex, selfenergy and box graphs with tb and g to the 
asymmetry A^ p at hadron level as functions of m H + are shown on Fig. 5. The 
large effect seen on the figure is mainly due to the phase of A t , and the asymmetry 
reaches its maximum for a maximal phase <pA t = vr/2. The phase of Af, does not 
have a big influence on the asymmetry and therefore we usually set it zero. The 
four kinks in all three lines denote the thresholds tuh+ = m>ii + m l >*> J = 1> 2, see 
again Table 1. 

The asymmetry Ap P reaches its maximum value at tan/5 = 5 and falls down 
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Figure 4: The contributions to the asymmetry Ap at parton level for the chosen 
set of parameters (59) as a function of y/1, m H + = 700 GeV. The red dotted line 
corresponds to box graphs with a gluino, the solid blue one to the vertex graph 
with a gluino, the green dashed one to the — if ± selfenergy graph with a tb 
loop, and the pink dash-dotted to the box and vertex graphs with 




Figure 5: The contributions to the asymmetry Ap P at hadron level for the chosen 
set of parameters (59) as a function of m H +. The red dotted line corresponds to 
box graphs with a gluino, the solid blue one to the vertex graph with a gluino, and 
the green dashed one to the W ± — H ± selfenergy graph with a tb loop. 

quickly with increasing tan /3. This dependence for m#+ = 550 GeV is shown on 
Fig. 6. 

In Fig. 7 we present the dependence of Ap P as a function of Mq(= = M^). 
The selfenergy contribution is first the biggest one, but it goes down to zero at 
Mq ~ 467 GeV because then the decay channel H + — ► closes. The kink at 
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Figure 6: The contributions to the asymmetry Ap P at hadron level for the chosen 
set of parameters (59), but m~ g = 450 GeV, as a function of tan /?, m#+ = 550 GeV. 
The red dotted line corresponds to box graphs with a gluino, the solid blue one to 
the vertex graph with a gluino, and the green dashed one to the W ± — H ± selfenergy 
graph with a tb loop. 

Mq ~ 450 GeV denotes the threshold H + — > t\b%. The box contribution has its 
maximum of ~ 12% at the threshold of H + — > t\bi 

A C P P 0.05 
0. 

-0.05 
-0.1 

350 400 450 500 550 600 

Mq [GeV] 

Figure 7: The contributions to the asymmetry Ap P at hadron level for the chosen 
set of parameters (59), but nig = 435 GeV, and m#+ = 800 GeV, as a function of 
Mq. The red dotted line corresponds to box graphs with a gluino, the solid blue 
one to the vertex graph with a gluino, and the green dashed one to the W ± — H ± 
selfenergy graph with a tb loop. 

In Fig. 8 the dependence of the three leading contributions to Ap^ p as a function 
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of nig is shown for m#+ = 550 GeV. Of course, the selfenergy contribution is 
independent of the gluino mass, being about —8 %. The vertex contribution has a 
maximum and the box contribution a minimum at m g ~ 425 GeV, and then their 
absolute values decrease. The box contribution is at the minimum about —27 %. 




_q 3 D l l l I I l l I l l l I I I l L 

' 400 600 800 1000 1200 
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Figure 8: The contributions to the asymmetry Ap P at hadron level for the chosen 
set of parameters (59) as a function of m g , m H + = 550GeV. The red dotted line 
corresponds to box graphs with a gluino, the solid blue one to the vertex graph 
with a gluino, and the green dashed one to the W — H ± selfenergy graph with a 
tb loop. 



5.2 Production and decay asymmetry 

First we want to add a few remarks on the branching ratios (BR) of the relevant 
decays. Fig. 9 shows the tree- level BRs of H + as functions of mu+- For small 
nin+, below the tb threshold, the dominant decay mode is H ± — > tb, with BR 
Ri 1, while the BR of — > vt^ is of the order of a few percent, decreasing with 
increasing . When the H ± — > tb channels are kinematically allowed, they start 
to dominate [21], and the BR of H ± — > z/r ± to a good approximation becomes zero. 
However, the BR of — > tb remains stable of the order of 15 — 20 %. The BR of 
— ► W ± h° reaches a few percent for small tan/3 in a relatively narrow range of 
nin+ [10]. In the considered range of parameters this decay is very much suppressed 
and we do not investigate it numerically. 

In Fig. 10 we show the total production and decay asymmetry A^ p at hadron 
level, for f = tb and / = z/r ± . Though for H ± — > z/r ± it can go up to ~ 20% for 
nin+ ~ 650 GeV, the BR of this decay in this range of H + masses is too small and 
observation at LHC is impossible. 
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Figure 9: The tree- level branching ratios of H + for the chosen set of parameters 
(59) as a function of m H +. 



Figure 10: The total asymmetry A CP at hadron level for the chosen set of param- 
eters (59) as a function of m#+. The blue line corresponds to the case where H ± 
decays to tb, and the green one to H ± decay to tu t . 

As we have shown analytically, the total asymmetry in the production and decay 
is approximately the algebraic sum of the asymmetry in the production A^, p , and 
the asymmetry in the decay A^f ' . One would think that the total CPV asymmetry 
will be large. Moreover, the CP-asymmetry in the decay alone is large [8]. 

Let us consider the subsequent decay H ± — > tb decay. In this case the contribu- 
tions coming from the selfenergy graph with tb in the loop in the production and 
in the decay exactly cancel. This cancellation occurs in general for any possible 
if ± selfenergy loop contribution to the s- or t-channel. It can be easily shown by 
writing down the matrix element for the whole three particle final state process. As 
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an illustrative example, let us consider the contribution of both selfenergy graphs 
with tb from the production and from the decay, to the s-channel of our process, 
see Fig. 11. The matrix element, representing the sum of the two graphs in Fig. 11 
reads 

2 | 



M 3, s = _Ws_l_ u{ph){yt p R + yb p L ) u{ - Pt ) [G,) l3 V} Ll Tl% + (G^KftRbj 

-« s (pt)7 PLp b T sr -f u r {p b )e x {p 



s 2 

. « + - "Wo* - »|)(p| - <") ^'^^^^ ■ 

(60) 

with J g = J d D q/(2ir) D , pf. = p H + + g and = g. It is clearly seen that (60) 
contains the sum of the couplings and their complex conjugate ones as a common 
factor. The expression in the box brackets at the end of the first row can be written 

as 

(GdijKiK) + (GdhKiKj = 2Re[(G 4 ) ij v} Li V} Lj } , (61) 

i.e. in this case the imaginary part of the couplings cancels. As the presence of a 
non zero imaginary part of the couplings is necessary for having CPV, it is clear 
that in this case the contribution of the selfenergy graphs on Fig. 11 to the CPV 
asymmetry is exactly zero. Our numerical study shows that the contributions of the 
vertex graphs from the production and from the decay also partially cancel with the 
box diagrams contribution. However, as the box graphs do not have an analogue 
in the decay, their contribution remains the leading one, see Fig. 12. 




Figure 11: The selfenergy contributions with tb in the loop to the s-channel of the 
considered process, from the production and from the decay. 

Fig. 13 shows the asymmetry in the production, Ap P , in the decay of H ± — > tb, 
A£) P tb , and the combined one, Af b p = Ap P + Ap) P tb , as a function of the phase 
(f>A t - All three curves are symmetric for (pA t —<fiA t because 0^ = 0. Ap P and 
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Figure 12: The contributions to the total asymmetry Af b p at hadron level for the 
chosen set of parameters (59) as a function of rriH+, with the same parameter set 
as used in Fig. 5. The dotted (red) line corresponds to box graphs with a gluino, 
the axially symmetric solid and dashed (green) lines correspond to the — H ± 
selfenergy graph with a tb loop in the production and the decay, respectively, and 
the other two solid and dashed (blue) lines correspond to the vertex graph with a 
gluino, again in the production and the decay, respectively. 



A^) P tb have negative relative signs with a maximium/minimum at \<f>A t \ ^ 0.457T, 
with \Ap P \ ~ 20% and |^4^ b | ~ 16% there. Due to this cancellation the resulting 
asymmetry A^ p is reduced to < 4%. 




Figure 13: The asymmetry A CP at hadron level for the production and the decay 
into tb only, and the total one, for the chosen set of parameters (59) as a function 
of (p At . 
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In order to avoid the cancellation, we can study the mass range of H ± , before 
the £i&i channel opens. In this case, having in mind our results in [8], the CP 
effects in the decay will be negligible and CPV will arise mainly in the production 
process due to the vertex and box contributions with g in the loops. For instance, 
for m H+ = 400 GeV and m g = 450 GeV we get Ag p = -3%. 

In Fig. 14 the total asymmetries Af b p and A„ p as functions of m g are shown for 
m H+ — 550 GeV. Both asymmetries are negative. They have their largest values 
at m g ~ 425 GeV, with A° p ~ -27 %, Af b p ~ -12 %, and then they decrease to 
zero. 




Figure 14: The total asymmetry A CP at hadron level for the chosen set of pa- 
rameters (59) as a function of m g , mn+ = 550 GeV. The (red) line denoted by tb 
corresponds to the case where if ± decays to tb, and the (green) one denoted by tv t 
to H ± decay to tu t . 

In Fig. 15 the dependence of Af b p on the absolute value of A t is shown for three 
different values of m H +. 

As already mentioned in the introduction, the phase of \i is strongly constrained 
by the measurements of the EDMs. Nevertheless, we also studied the dependence 
of A CP on 0^. Using /i — 700 e l z GeV instead of /i — — 700 GeV, we get for 
m- = 450 GeV in Fig. 14 the asymmetry Af b p ~ -22%. 

The tan/3 dependence is shown on Fig. 16, for f = tb and / = tv t . In both 
considered cases the asymmetry A CP has its maximum at tan/3 ~ 5, with A^ p ~ 
12% and A„ p ~ 21%. It approximately vanishes for tan/3 > 15. 

We have compared our results with those in [13, 14]. Our numerical results are 
in good agreement with [13], but we disagree analytically and numerically with [14], 
where in addition the box contributions are missing. 

The production rate of H + at the LHC for m#+ = 550 GeV and tan/? = 5 is 
~ 15 fb (including a K-factor from QCD of 1.5, see [22]), and the BR(H + tb) 
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Figure 15: The total asymmetry Af b p at hadron level for the chosen set of parame- 
ters (59) as a function of \A t \, for three values of (in GeV). 
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Figure 16: The total asymmetry A CP at hadron level for the chosen set of param- 
eters (59) as a function of tan (3, m g = 450 GeV, m H + = 550 GeV. The (red) line 
denoted by tb corresponds to the case where if ± decays to tb, and the (green) one 
denoted by tu t corresponds to H ± decay to tu t . 

~ 20%. For an integrated luminosity of 200 fb" 1 we get N ~ 600 and = 25. For 
rrig = 450 we get Af b p ~ 0.12 and therefore the statistical significance ^/N Af b p ~ 4. 
But because of the large background, the actual signal production rate will be 
most likely reduced and the statistical significance might be too low for a clear 
observation in H^t production in the first stage of LHC However, at SLHC with a 
design luminosity bigger by a factor of ~ 10, such a measurement would be worth 
of being performed. 
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6 Conclusions 



The MSSM with complex parameters in particular with A t complex, gives rise to 
CP violation in the production of H ± , pp — * H ± t + X, and in the decays of H ± 
to tb, v T r ± and W^h at one-loop level. We have calculated the corresponding 
asymmetries between the H + and H rates both in the production and in the 
decays. A few improvements have been made with respect to previous calculations. 

We have performed a detailed numerical analysis studying the dependence on 
the important parameters. A peculiarity is that in the case of pp H ± t + X with 
if ± — > tb, the contributions coming from the selfenergy graph with tb in the loop in 
the production and in the decay exactly cancel. Nevertheless, the asymmetry can 
go up to ~ 12% for m H + ~ 550 GeV, mainly due to the box graphs with gluino in 
the production process. A measurement of the CP violating asymmetries at LHC 
can give important information of the parameters of the MSSM, especially on A t 
and its phase. 

Acknowledgments 

The authors acknowledge support from EU under the MRTN-CT-2006-035505 net- 
work programme. This work is supported by the "Fonds zur Forderung der wis- 
senschaftlichen Forschung" of Austria, project No. P18959-N16. The work of E. C. 
and E. G. is partially supported by the Bulgarian National Science Foundation, 
grant 288/2008. 

A Masses and mixing matrices 

The mass matrix of the stops in the basis (t^, Ir) reads 




We have 




(63) 



Analogously, the mass matrix of the sbottoms in the basis (&£, 6r) 

2 / M? - m| cos 2f3(\ - \ sin 2 6w) + m l (^b ~ A* tan P) m b 

^ l ~ \ (A b - /j,* tan (3) m b - \m\ cos 2/3 sin 2 Ow + 



.(64) 
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is diagonalized by the rotation matrix lZ b such that TZ b ^Ai^TZ b = diag(m~ , m~ ). 

ft 6 has the same structure like 1Z 1 given with (63), and one can obtain it by making 
the interchange t — > b. 

B Interaction Lagrangian 

The part of the MSSM interaction Lagrangian used in our analytical calculations 
is given in this section. 

The interaction of the charged Higgs boson with two quarks reads 

C H±tb = H + 1 (y t P L + y b P R ) b + H-b (y b P L + y t P R ) t , (65) 

where the P L and P R are the left/ right projection operators 

^ = ^(1-75), P fl = ^(l + 76), (66) 
y t and y b are the tree-level couplings 

y t = ht cos f3, y b = h b sinP, (67) 
with h t and h b - the top and bottom Yukawa couplings 

ht — — 7= 7; 5 hb = —}= 77 • (68) 

\J1m w sin jj \J1m w cos p 

The interaction of two quarks with gluon exchange is given by 

C qqg = -g s T>> G; q t Y q s , (69) 

with i,j = 1, 2, 3 and a = 1, .., 8. 

The interaction of the charged Higgs boson with two squarks is described by 

Z&uh = ( G *h H+ 1 * h + H ~ ~ b *> *i . ( 7 °) 

with i,j = 1, 2, 

G A = n^G A n\ (71) 

and the matrix G4 is given by 

(h b m b sin (5 + htirit cos (3 — \/2 g mw sin /3 cos j3 h b (At sin j3 + fi cos (3) \ 
■ (72) 
h t {A t cos (3 + fi* sin f3) h t m b cos (3 + hbirit sin /? y 

The W^ ± -squark-squark interaction Lagrangian reads 
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where i,j = 1, 2 and A B = A {d„B) - (d^A) B. 
The squark-quark-gluino interaction is given by 

C q¥g = -V2g s T° k [UKl\e-^P L -K%e^P R ) q k qf 

+ q^TZl ei**P R - K% e-^P L ) ~g a $ } , (74) 

with i — 1, 2, j, k — 1, 2, 3, and a = 1, ..8. 

The interaction of two quarks with W-boson exchange is described by 

Cw±tb = -^( W ^^ P Lb + W~b^P L t). (75) 

C Passarino-Veltman integrals 

The definitions of the Passarino-Veltman two-, and three-point functions [23] in 
the convention of [24] and the derived analytical expressions for their imaginary 
parts [10, 25] are given in this section. 

The PV two-point functions are defined through the 4-dimensional integrals, as 

Boiplmlml) = l^JdPq^p, (76) 

B^p\,ml,ml) = Jd D q^^ I =p lt ,B 1 {p 2 1 ,m 2 ,m 2 1 ), (77) 

where we use the notation 

V° = q 2 - m 2 and V j = (q + Pj ) 2 - m) . (78) 

In the rest frame system of the (decaying) particle with impulse p 1 , for the 
imaginary part of B we get 

im^.mj.^-'^y^, (79) 

where the A-function is defined as 

X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz. (80) 
In order to derive the imaginary part of B\ we use the relation 

2k 2 B x {k 2 , m 2 , ml) = A (m 2 ) - A (m\) + [m\ - m 2 - k 2 )B (k 2 , m 2 , m\) . (81) 
Having in mind that lmA (m 2 ) = 0, we obtain 

Tm R (M 2 m 2 m 2 \ - 7r ( m l ~ m ~ Mf)\ 1 / 2 (Mf, m„, m\) 
VmB 1 {M 1 ,m G ,m 1 ) . (82) 
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The PV three-point functions are defined as 



1 



ITT" 



d D q 



£>0£>1£>2 ' 



cM, (pi - P2) 2 ,P2> m o> m l) = A / ^ v^v 2 = Plf * Cl + P2f * C2 ' 

(83) 

For the absorptive parts of the integrals C , C\ and C 2 in (83) in the rest frame of 



Pl~P2,M 3 




q,m 



P2,M 2 



Figure 17: The selfenergy and vertex type diagrams corresponding to the PV two- 
and three-point integrals, when particles with masses m and m\ are on mass shell. 

the (decaying) particle with momentum pi, when particles with masses m and mi 
are on mass shell, Fig. 17, we obtain 



hnC (MlMlMl,m*,mlml) = 
Im C^M 2 , Ml Ml ml m 2 1: m 2 ,) = 

hnC 2 (MlMlMl,m*,mlml) = 
Here we have 

A = M\Ml - ( Pl p 2 ) 2 , 



7T 



AV2( Ml 2 ,M|,M|) 

M 2 A-( Pl p 2 )B 
A 

( Pl p 2 )A-M 2 B 



In 



,a + b. 
a — b 



a = Ml + m„ - m\ + 2 g ^ , 



M 2 + M| - Mf 
2 



Tx{m\-m\- M\) a + b 
~ ~2\ 1 / 2 {M 2 ,M 2 ,M 2 ) n| ^6'' 
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7T 



A 1 /2(M 1 2 ,M 2 2 ,M|) \2 



J(M| + m„ - ml) In' £ ' 



a — b 



+ 2 | g | |p 2 



_ (m 2 - m 2 - M 2 )(M 2 + M 2 - Mf) 

H F2 ~ 



9 I |P2 



4M 2 

A 1 / 2 (M 2 , m 2 , m 2 ) A 1 / 2 (M 2 , M 2 , Mf) 
4M 2 



(84) 

(85) 
(86) 



(87) 
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In (84)-(86) and further, we denote masses of the on-shell particles with bold font. 
Generally, the functions C , C 1 and C 3 have three terms according to the three 
possible cuts over the particles in the loops, and which of them is non zero depends 
on the kinematics. The imaginary part of Co reads 

Im C (M 2 , Ml Ml m 2 , m\, m\) = Im C (M 2 , Mf , M 2 , mg, m 2 , m 2 ) + 

+Im C (M 2 , Ml Ml ml, m 2 , ml) + Im C (M 2 , M 3 2 , M 2 2 , mg, m 2 , ml) . (88) 

In terms of the analytic formulas (84)-(86) we write 

Im C (M 2 , Ml Ml ml, m\, m\) = Im C (M 2 , M 3 2 , M 2 , mg, m 2 , m 2 ) + 

+Im C (Ml Ml Ml m 2 , m 2 , mg) + Im C (M 2 , M 3 2 , M 2 , mg, m 2 , m 2 ) . (89) 

Im d(M 2 , M 3 2 , Ml ml, m\, m 2 ) = Im d(M 2 , M 3 2 , M 2 2 , mg, m?, m 2 ) + 

+Im d(M 3 2 , M 2 , M 2 , m 2 , m?, m 2 ) + Im C 2 (M|, M 3 2 , M 2 , mg, m 2 , m 2 ) , (90) 

Im C 2 (M 2 , M 3 2 , M 2 , m 2 , m\, m\) = Im C 2 (M 2 , M 3 2 , M 2 , m 2 , m?, m 2 ) + 
+ [ImC (M 2 , Ml Ml ml, m 2 , m 2 ) + Imd(M 3 2 , M 2 , M 2 2 , m 2 , m 2 , m 2 ) + 
+Im C 2 (Ml Ml Ml m 2 2 , m?, m 2 )] + Im d(M 2 2 , M 2 , M 2 , m 2 , m 2 , m 2 ) . (91) 
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